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ABSTRACT 

The drag force on a satellite of mass M moving with speed V in the gravitational field of a spherically 
symmetric background of stars is computed. During the encounter, the stars are subject to a time-dependent 
force that alters their equilibrium. The resulting distortion in the stellar density field acts back to produce 
a force Fa that decelerates the satellite. This force is computed using a perturbative technique known as 
linear response theory. 

In this paper, we extend the formalism of linear response to derive the correct expression for the 
back-reaction force Fa that applies when the stellar system is described by an equilibrium one-particle 
distribution function. Fa is expressed in terms of a suitable correlation function coupling the satellite 
dynamics to the unperturbed dynamics of the stars. At time t the force depends upon the whole history of 
the composite system. In the formalism, we account for the shift of the stellar center of mass resulting from 
linear momentum conservation. The self-gravity of the response is neglected since it contributes to a higher 
order in the perturbation. Linear response theory applies also to the case of a satellite orbiting outside the 
spherical galaxy. 

The case of a satellite moving on a straight-line, at high speed relative to the stellar dispersion velocity, 
is explored. We find that the satellite during its passage rises (a) global tides in the stellar distribution and 
(b) a wake, i.e., an overdense region behind its trail. 

If the satellite motion is external to the galaxy, it suffers a dissipative force which is not exclusively 
acting along V but acquires a component along R, the position vector relative to the center of the spherical 
galaxy. We derive the analytical expression of the force, in the impulse approximation. 

In penetrating shortlived encounters, the satellite moves across the stellar distribution and the wake 
excited in the density field is responsible for most of the deceleration. We find that dynamical friction rises 
from a memory effect involving only those stars perturbed along the path. The force can be written in terms 
of an effective Coulomb logarithm which now depends on the dynamical history, and in turn upon time t. 
In A is computed for two simple equilibrium density distributions: It is shown that the drag increases as the 
the satellite approaches the denser regions of the stellar distribution. The braking force then stays almost 
constant. When the satellite crosses the edge of the galaxy the force does not vanish but declines with time 
as R ^ oo. 

In the case of a homogeneous cloud, we compute the total energy loss. In evaluating the contribution 
resulting from friction, we derive self-consistently the maximum impact parameter which is found equal to 
the length travelled by the satellite within the system. Tides excited by the satellite in the galaxy reduce 
the value of the energy loss by friction; in close encounters this value is decreased by a factor ~ 1.5. 
Subject headings: galaxies: clustering - stars:stellar dynamics 



I. INTRODUCTION 



A massive object M moving with velocity V through a background of field stars suffers a dissipative 
force known as dynamical friction. In the original formulation of Chandrasekhar (1943), this force is the 
result of the momentum exchange between the particle M and the stars of the background. In a uniform 
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isotropic stellar system, the uncorrelated superposition of these binary encounters leads to a frictional force 



where p is the mass density of the stars having speed less than and hmax and hmin are the maximum and 
minimum impact parameters for which encounters can be considered effective, hmax is conventionally set 
equal to the size of the stellar system surrounding the object (to avoid a divergence) and bmin is the larger 
between the characteristic radius of the satellite M and the impact parameter for a 90° scattering event 
(White 1976). In the derivation of the above equation the stars interact solely with the gravitational field of 
the incoming particle and move along trajectories that arc straight lines at infinity. Because of the isotropy 
of the surrounding stellar field the force on M, responsible for its deceleration, is acting exclusively along 
the direction of motion, i.e., along V. Dynamical friction can be viewed, alternatively, as a force originating 
from the overdensity excited in the stellar field by the particle in its motion: a wake develops behind M 
which is causing its deceleration (Mulder 1983). 

Dynamical friction plays a central role in many astrophysical environments: It intervenes, for example, 
in the evolution of galaxies accreting small dwarfs (Lin & Tremaine 1983; Quinn & Goodman 1986; Binncy 
& Tremaine 1987), in galaxy cannibalism, in the collision and merger of galaxies (Barnes & Hcrnquist 1992), 
and in the inspiral of a massive black hole at the center of dense stellar clusters (Bcgclman, Blandford & 
Rees 1980). The expression of the frictional force in a homogenous background has therefore been used to 
derive simple analytical estimates of the orbital decay times and to track evolutionary scenarios, for these 
cases. In these systems however, the ambient medium is neither uniform, nor infinite and this calls for a 
deep analysis of this process. 

Since the work by Chandrasekhar, dynamical friction has been studied using various techniques. Tremai- 
ne and Weinberg (1984), and Weinberg (1986) studied the loss of angular momentum suffered by a satellite 
moving in the gravitational field of a spherical galaxy, along a circular orbit. Using a perturbative approach 
(that refers also to works by Lynden-Bell & Kalnajs 1972) for the description of the interaction of the satellite 
with the background of self-gravitating point masses, they found that angular momentum is exchanged 
secularly with those stars whose orbits are commensurate, i.e., in near-resonance with the satellite's orbit. 
In applying the technique to the case of a satellite (modelled as a Plummer sphere) orbiting inside of a singular 
isothermal sphere, Weinberg (1986) found that the resulting torque on M is equivalent in magnitude to the 
one obtained using Chandrasekhar formula, applied to a locally homogeneous system: The maximum impact 
parameter hmax is found to be comparable to the size of the satellite orbit Rs- This result suggests that the 
braking of a satellite depends more on the local properties of the stellar background, despite the long-range 
nature of the gravitational force. 

For a parent galaxy modelled as a polytrope and a satellite as a non-deformablc Plummer sphere similar 
findings were inferred by Bontekoe & Van Albada (1987). In a self-consistent A'^-body numerical calculation 
these authors observed that the decay process can still be described as if it were a local phenomenon in which 
the density field around the object determines primarily the magnitude of the torque and in turn the decay 
time, for a satellite moving along a circular orbit (see also Farouki, Duncan & Shapiro 1983). Bontekoe 
& Van Albada (1987) also explored the evolution of a satellite moving around the primary galaxy on a 
grazing orbit, and found that after a few revolutions the satellite, losing progressively energy and angular 
momentum, enters the stellar distribution and sufi^ers complete merger. This case illustrates the inadequacy 
of Chandrasekhar formula since it would predict no drag for a satellite orbiting outside the outer edges of 
the companion galaxy. But not only, it clearly indicates that the satellite experiences a deceleration that 
rises in response to the perturbation that is excited in the stellar field during the interaction. This is a global 
phenomenon that can be explained in terms of an interaction with resonant stars in the galaxy (Tremaine 
& Weinberg 1984). Thus a continuity should exist in the physical processes inducing orbital decay. 

The aim of this paper and of the accompanying one is to clarify the role of the tidal and frictional 
forces in affecting the motion of the satellite. This is achieved within the context of the theory of linear 
response. The method is close to that explored by Bekenstein and Maoz (1992) who relate the frictional drag 
to the stochastic components of the background forces, thus establishing a connection between dynamical 
friction and the fluctuation-dissipation theorem. Their method proved useful in giving a unified picture of 
the process of energy loss and heating of a test particle traveling through a uniform stellar system. In a 
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subsequent work, Maoz (1993) provided an approximate expression for the energy loss experienced by a 
massive particle moving in a nonhomogeneous medium. He noted that, owing to the nonuniformities present 
in the underlying stellar system, dynamical friction no longer depends on local background characteristics 
and that the force acquires a component in the direction of the spatial density gradients. Nevertheless he 
does not compute the force nor give a complete and consistent description of the process. 

In this paper we compute the force of feed-back on a satellite moving in a background which is nonuniform 
and spherically symmetric, improving the perturbative technique first developed by Kandrup (1981,1983; 
see also the recent work by Nelson & Tremaine 1997). The coUisionless system of stars is described by 
the equilibrium distribution function /q. During the passage of the satellite M, the stars are subject to a 
time-dependent gravitational force and on account of that, the distribution function modifies: Its new value 
/ = /o + A/ is calculated using the theory of the linear response. The perturbation induced by the satellite 
acts back producing a force Fa which is computed, accordingly, as expectation value on the response function 
A/ of the microphysical force resulting from the interaction of M with the stars of the background. The 
characteristics of this force are closely related to the way the stellar background responds to the perturbation 
and depend, as in memory effect, on the whole history of the particle and stars since the moment in which 
the interaction is turned on. The formalism also applies to the case of a satellite orbiting outside of the 
stellar distribution. 

In this paper we wish to shed light into the nature of the force and derive a natural decomposition in 
terms of a global and local response of the background. This work is also preparatory to the analysis that 
will be carried in paper II, on the role of tides in causing the orbital decay of a satellite accreting onto a 
companion galaxy, a process relevant to the clustering of structures in cosmology. 

In addressing all these problems, for ease of analysis, and consistently with the hypothesis that stars in 
equilibrium can be considered as a coUisionless system, the equilibrium distribution function /o is factorized. 
This assumption is customarily used to compute mean properties of the response, as in Weinberg approach 
(1986). We here find that the introduction of such hypothesis cancels out artificially an important correlation 
present among the stars: Under the action of their self-gravity, the stars perform bound orbits around the 
barycenter of the galaxy. The displacement that the center of mass suffers during the encounter with the 
satellite (owing to linear momentum conservation) represents a collective response neglected in previous 
analysis. This response needs to be included for a correct estimate of the transfer of orbital energy into the 
internal degrees of freedom of the galaxy. 

In this paper we modify the formalism of the theory of linear response to account for this important 
effect. This modification explains the longstanding problem of the differences in the sinking times found in 
N-body simulations (White 1983; Zaritsky & White 1988; Hernquist & Weinberg 1988) and in semi-restricted 
methods (Lin & Tremaine 1983). We present here the method used to include this effect but defer to paper 
II for a more thorough analysis. 

The layout of the paper is as follows: In §2 we sketch the method used to compute the macroscopic force 
on M, using a derivation which is independent of that given by Kandrup (1981). In §3 we test the theory 
against known results, focusing on shortlived distant encounters. In §4 we derive the formal expression of 
the force Fa on M that accounts for the shift of the center of mass of the galaxy. We then compute in §5, 
the force on the satellite M moving at high speed through a nonuniform stellar background characterized 
at equilibrium by a maxwellian distribution function. In §6 we compute the frictional force on a satellite 
moving across a finite size cloud with given density profile, and compute the shape of the tidal distortion. 
In §7 we evaluate the extent of the total energy loss. In §8 we present our conclusions. 



II. THEORY OF LINEAR RESPONSE 



2.\Response Function and Back-Reaction Force 

Consider a particle (i.e., a satellite) of mass M and velocity V moving through a coUisionless system 
of N stars. The stars have mass m <^ M and are described by an equilibrium distribution function /q. At 
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R(t), the massive particle experiences an A/'-body force resulting from the gravitational interaction with the 
background stars 

N 



F{t) = -GMmY ^St^ 



(1) 



where rj is the position of the i— th star at time t. By action-reaction, each of the N stars of the background 
experiences a perturbing force 

Under the influence of the time dependent force Fj, the stellar distribution function / is altered and the 
distortion with respect to equilibrium can be described, to first order in the perturbation, in terms of the 
response function A/ such that 

f = fo + Af. (3) 

Because of the change in the distribution function induced by the satellite M, the force on M changes. This 
macroscopic force can be determined as mean on the distribution function / and is decomposed into two 
terms 

< F >/„+A/=< F >/„ + < F >A/ . (4) 

The force < F > originates from the interaction of M with the sm,ooth potential generated by the unper- 
turbed stellar system; the term < F >a/ instead will represent the force of back-reaction whose magnitude 
depends on the response function A/. As a first step, we calculate the perturbation A/ on the stellar 
distribution function. 

Let — F = J2iLi Fi be the applied time dependent force acting on the stellar background. The distri- 
bution function / = /g + A/ can be computed, to first order in the perturbation, using the theory of linear 
response (TLR hereafter). For this purpose we need to specify the form of the disturbance and the properties 
of the stellar system in the unperturbed state through the equilibrium Hamiltonian Hq. 

In the coUisionless stellar system 

N 2 
i=l 

where describes the smooth self-interaction, i.e. the mean field of the gravitating system. In presence of 
the time dependent external force induced by the satellite, the Hamiltonian modifies and the perturbation 
Hi is simply 

Hi results from the superposition of the gravitational force that the incoming body exerts on each particle 
of the system. Hi is regarded as an explicit function of t depending on the instantaneous position R(t) of 
the particle M and as a function of the 6N coordinates (r, p) of the phase space: hence, Hi = Hi{t; F). (F 
denotes shortly the phase-space coordinates) 

According to the Liouville equation, the evolution of the distribution function / is given by the equation 

^ = {Ho + Hi{t),m}. (7) 

If the distribution function /g depends on the phase-space coordinate only through Hq the Poisson bracket 
{Ho,fo} vanishes identically, and to first order in the perturbation, we can approximate equation (7) as 

^^M={Ho,Af} + {Hi{t),fo}. (8) 

The formal solution is 

Af{t;T)= f dse-'^'-^^"- {Hi{s),h], (9) 

J — OO 
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where C = i{Ho, } denotes the LiouviUe operator acting on the 6N coordinates of the phase space. In 
equation (9) the expression of the Poisson bracket reads 



(10) 



In defining the hmits of integration of equation (9), the external force is applied in the infinite past when 
the system was in a stationary state and the particle M far away: accordingly, Hi vanishes for f. — oo. 

Having calculated the distortion A/, (to first order in the perturbation Hi) we can compute the expec- 
tation value of the microphysical force F acting on the satellite M at time t : 



<F{t) >Af= J dT Af F{t;T). 
For a perturbation of the form given by equation (10) 

f dsfdVe--^- ^ ^^^^ 

,; _ 1 J —OO J 



< F(i) >A/= 



-GMi 



Vpjo' 



R(s)-ri|3 



e'*^ [F(t; r)] 



(11) 



(12) 



where we used the property that £ is a self-adjoint operator. In equation (12) the LiouviUe operator acts 
on the phase-space coordinates (F): applied to a generic dynamical variable A (representing any component 
of the phase-space coordinates which satisfies equation dA/dt = iCA ) C carries A along its evolution from 
time i = to time t , i.e., 



A{V{t)) = e'*^A(r(0)) 



(13) 



Hence, we interpret e^'^'^F at an arbitrary time t as the value of the force F acting on M at R(t) resulting 
from the gravitational interaction with the particles that evolved from t = to time t along a path whose 
dynamics is uniquely determined by the unperturbed Hamiltonian Hq. Likewise, e""^ acts on the scalar 
product to yield accordingly 



< F(i) >A/= -GMm 



„_1 J —OO J 



dT 



|R(s)-r,(s)|3, 



F{t;r{t)) 



(14) 



We will shortly denote < F{t) >a/= Fa hereafter. This equation provides the formal expression of the 
force acting on M in terms of the microscopic force F and of the equilibrium distribution function /q. We 
notice that within TLR, the response function A/ is of order "G" ; Hence, the resulting force describing the 
feed back mechanism acting on M is of order " " . Equation (14) contains all details of the microphysical 
interaction of the A'' stars with the perturber: It accounts for the self-gravity of the stellar equilibrium and 
for the actual dynamics of the perturber. 

The diSiculty in evaluating the multiple integrals of equation (14), involving the actual dynamics of the 
underlying stars, is enormous particularly if we want to use analytic tools to gain insight on the mechanism 
producing friction. To overcome this problem we introduce a major simplification knowing that the stellar 
motion in a galaxy can be described in terms of a mean field potential (Binney & Tremaine 1987): The system 
can therefore be regarded as collisionless and stars behave as independent particles. The corresponding 
distribution function /o(r) can therefore be written in terms of the one-particle phase space density /°P: 



fo = lllj°P{r,p). 



(15) 



Under this hypothesis, due to the statistical independence of the particles, the cross-correlation terms 
(i.e., those involving different indices in the summation of eq.[14]), cancel identically, in the limit N ^ 1. 
Only the term of self-correlation survives giving rise to a macroscopic force 



Fa = [GMfNm'^ f ds f dardap 

J — OO J 



Vp(«) f 



op 



R(s) - r(s) 



|R(s)-r(s)|3 



\m-r{tw 



(16a) 
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where the phase-space vohime d^r d^p and /°p arc independent of time and can be referred to time t. An 
equivalent expression for the force can be obtained by applying the evolution operator g-^i*-^)^ to F{t; F) 
in equation (12) and reads 



Fa = [GM]27Vm2 / ds f rfsrrfaP 

•/ — OO •/ 



|R(a) 



R(i) - r(i - s) 
|R(i)-r(t-5)|3 



(166) 



(equation (16b) will be used in paper II). 

Introducing the velocity vector v = p/m, we can express /°p as a function of (r, v). If the equilibrium 
system is described by a Gaussian distribution function 

\ 3/2 1 



/°P(r,v) = no(r) j exp(--/3mi;^) (17) 

we have Vv = —/3,f°^ v. In equation (17), r and v denote the absolute values of r and v and the coefficient 
P = (m(7^)~^ is directly related to the one-dimensional dispersion velocity a characteristic of the equilibrium 
stellar distribution; no{r) is the unperturbed probability density which for a spherically symmetric system 
is a function of r only. /°p is defined so that / dar no(r) = 1. With the change of variables (r.p) (r, v) 
and according to the normalization used we recall that (i3r(i3p/°P(r,p) is equal to dsr (i3v/°P(r, v). 

The force Fa at time t is found to depend on the whole " history " of the composite system, i.e., on 
the dynamics of the satellite, and on the dynamics of stars of the background as determined by the (mean 
field) Hamiltonian Hq. If we introduce the tensor 

T = -[GM] Nm /3|R(3)_r(,)|3|R(i)_r(t)|3 (l^) 

correlating the microscopical force at time s with that at time t, we can interpret equation (16) written in 
the form ^ 

Fa = J ds J rfsr dgv f°Pv''{s)T'"' (19) 

as a manifestation on the fluctuation-dissipation theorem relating the averaged force of back-reaction to the 
time integral of a suitable dynamical correlation function of the unperturbed system (Bekenstein & Maoz 
1992; Nelson & Tremaine 1997). Due to the long range nature of the gravitational interaction all stars in 
the galaxy give a sizable contribution to the force. Fa is a force which is intrinsically nonlocal both in 
phase-space and time. In equation (19), v^T^'^ are shorthand for the summations (hereafter). 

Notice that in the steps that lead to equation (19) the perturber's motion needs not necessarely to be 
confined within the stellar distribution. We can therefore evaluate, within TLR, the force on the massive 
particle M when it is orbiting outside the galaxy. 



2.2 Force and Density Distortion 



Fa can be alternatively viewed as a force rising in response to the distortion caused by the satellite on 

the stellar density field. The satellite M induces a disturbance in the density distribution function n(r) that 
can be readily evaluated, within TLR, from equation (9). Relative to its unperturbed value no, the density 
changes by an amount equal to 

An(r) = j dar" d^v S^ir" - r)A/(r", v) (20) 
where 5^ denotes the 3-dimensional Dirac function. From the solution (9) and equation (10) we have 

An(r) = -GM £^ ds J dsv" dsyrSs{r"it - s) - r)V.nr", v) ■ ^^^^^^ ~ (^1) 
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In the stops that load to equation (21) wo havo usod oxphoitly tho solf-adjointnoss proporty of tho LiouviUe 
operator and the factorization of /o { C is apphed to the phase-space coordinates of the Dirac function). 
Under the disturbance (21), the satellite experiences a force 

Fa = -GMNm j dsvAnir)-^^^ (22) 

which coincides with equation (16b). The expression of the force Fa will be modified in §4 to account 
properly for the shift of the centre of mass of the stellar distribution, that occurs during the gravitational 
interaction. 

Given Fa, we can formally compute the total energy loss suffered by the satellite 

/ + 00 1 /■ 

dt Fa-V = --[GMmfNp j dsrdav (23) 

where 

.+00 R-r 

Hereafter R, r, and v will denote the vectors at the current time t and the volume element d^rd^v can be 
referred only for convenience to time t. (The derivation of equation (23) is summarized in Appendix A.) 



III. TESTING THE THEORY 



Kandrup (1983) vorifiod, in his oarly work, that equation (16) roproducos tho expression of the drag 
force derived by Chandrasc^khar for the case of a point mass moving through an infinite homogeneous stellar 
background. This important limit is recovered within TLR in the hypothesis that the stars and the perturber 
move along straight lines (we do not sketch here the derivation of this limit but refer to Kandrup). According 
to the original formulation (Chandrasekhar 1943), the gravitational interaction among the stars was neglected 
and the braking force on M was derived considering the momentum exchange resulting from the incoherent 
superposition of the binary encounters with single stars moving along hyperbolic orbits relative to M. TLR 
reproduces the correct results to order G^, giving in the high speed limit, i.e, for V ^ a a, frictional force 

Fdf = -4:Tt ' ' Nmnoln — (25) 



\bminj V 

where bmax is set equal to the characteristic length of the stellar system to avoid a "infrared " divergence. 
In TLR, the above expression is derived introducing also at short distances a "ultraviolet" cutoff, bmm, the 
minimum impact parameter for collisions at large scattering angle. The presence of a divergence at small 
relative distances between M and the stars, in our formulation, is a consequence of the approach used. 
Since we have treated Hi as a small perturbation, equation (9) becomes invalid when the distortion in the 
background produces a wake with overdensity An/no ~ 1. The density enhancement An/no (computed 
using equation (9) and neglecting the self-gravity of stars) 

- = ^m^ei»'-^'Herta+l] (26) 
no cr^ |R— r| 

(with a = (V/-\/2c7) • (R — r)/|R — r|, and x = Vl\/2a) becomes close to unity when the relative distance 
between a star and the satellite M is 



Max(a2,y2) 
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One can recognize that at such a relative distance the mean squared velocity of a typical field particle 
would increase by an amount of order a (V, in the high speed limit). For an extended (nondeformable) 
satellite e is instead set equal to physical extension of M. Thus e gives the "permitted" size the satellite 
(White 1976). In the expression of the force Fa the domain of integration is thus limited to |R — r| 3> e. 

Can the theory of LR be tested against other known results ? The energy loss of an object of mass M 
moving at high speed V relative to a galaxy characterized by a mean-square radius < > is known to be 

This equation is derived in the hypothesis that the distance of closest approach h largely exceeds the radius 
of the disturbed galaxy and that the effective duration of the encounter Tg be shorter than the crossing 
(dynamical) time of the stars Tdyn ~ (< > / < t;^ >)^/^ where < > is the mean squared velocity 
(< >= 3(7^). This limit is known as impulse approximation and applies when the perturber moves at 
high speed {V ^ a) relative to the stars. The energy loss (28) is of order "G^"and was derived by Spitzer 
(1958) on the basis of simple arguments, without direct knowledge of the force causing such a loss (Binney 
& tremaine 1987). 

Using equation (23) we recover Spitzer's results as in high speed encounters the energy loss of the satellite 
is exceedingly small (scaling with V^^ ), and TLR applies. The motion of the perturber can be approximated 
as linear while stars barc;ly move since the interaction is shortlived. Under these simplifications, equation 
(23) for the energy loss of M can be evaluated analytically (the calculation is sketched in Appendix A) to 
give 



= -2^- W - < > W + « [-1^ ) (29) 

AE as expected is the sum of two contributions. The first term, according to linear momentum conservation, 
is due to recoil, i.e., the effect representing the exchange of energy between M and the galaxy viewed as 
pointlike. The second instead describes the transfer of orbital energy into the internal energy of the stellar 
distribution. 



IV. TLR REVISED: THE FORCE IN THE GALAXY REST FRAME 



4.1 Correlation and shift of stellar bary center 



In this paper, we wish to gain insight into the mechanisms that ultimately cause the transfer of orbital 
energy into the internal degrees of freedom of the stellar distribution. For this purpose, we need to isolate 
from the force Fa those terms responsible to the recoil, representing only a global displacement of the 
collection of stars in the system with no relation to the "binding" mechanism that M experiences as a result 
of friction. We note that equation (16) derived under the hypothesis that stars, in equilibrium, behave as 
independent partieles, does not correctly describe the energy transfer during the interaction. At the basis of 
this inconsistency is the simplifying assumption of factorization of /, the equilibrium distribution function 
(eq. [15]). According to this hypothesis, the dynamics of each star is uncorrelated from the dynamics of 
any other star in the galaxy. In this approximation each star moves in the averaged potential generated 
by all the other stars. For a spherical distribution this potential gives rise to a central force binding the 
stars to the center of mass of the galaxy. In the encounter with a satellite, the center of symmetry of 
the stellar distribution suffers a rigid displacement due to linear momentum conservation. The stars thus 
coherently move responding to the satellite perturbation with a coordinated shift of their orbits. However, 
this correlation (present in the full formalism leading to eq. [14]) is lost, since the mean field dynamics we 
use to represent the unperturbed motion of the stars (eq. [15]) does not allow for the shift of the center of 
mass; the dynamics has lost its translational invariance. 

As a consequence, equation (23) for AE gives a consistent estimate of the energy loss only in three 
limiting cases; when (i) the satellite moves through a homogeneous infinite background of non-interacting 
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stars, (ii) when the interaetion is shorthved, as in the impulse approximation or (iii) when the center of 
symmetry of the galaxy is pinned. In the encounter of a satellite with a spherical galaxy this last hypothesis 
is clearly unphysical. The above considerations rise therefore a problem of consistency which can not be 
overlooked. 

To solve for this difficulty, we examine the process of dynamical friction not from the laboratory system 
but from the noninertial frame of reference comoving with the center of mass of the perturbed galaxy: 
X = (l/N) Newton's equations, in this frame, read 

M^^^^^=<F>/op+Fa (30) 



m^^^^^ = ^internal + (-F^) - ^ ^(-F,). (31) 

The satellite simply acquires a new mass equal to the reduced mass of the system /x = MNm/{M + Nm). 
Instead, the rhs of equation (31) for the motion of the i— th star acquires a new term, representing the 
opposite of the "mean" force exerted by the satellite on the galaxy, as a whole (Fj is the force acting on M 
due to the i— th star; eq. [2]). 

Guided by equations (30) and (31), we are led to modify Fa introducing a new correlation tensor T 
obtained subtracting such a " mean" force to the kernel [R — r]/|R — r|^ along the complete history until 
time t, yielding 



Fa = - [GMfNrn^li f ds f rfardsv/^P 

J —OO J 



v(s) 



R(s) - r(s) 
|R(s) - r(s)|3 



ci3r'no(r') 



R(,s) 



R(s) 



r-'|3 



R 



IR- 



IS- 



(32) 



Fa represents here the response force acting of M at time t as measured in the frame comoving with the 
perturbed stellar distribution. We now regard R and r as coordinates relative to the center of mass of the 
perturbed stellar system. 

In this frame, the disturbance in the density field is calculated according to equation (32), by introducing 
the modified kernel. For a Gaussian distribution function, we have 



An{r) =GMm(i (^^^ ' J ds j d3r" d3V S3{r"{t - s) - r)no{r")e-^i^"'''"^ 



R(g) - /• , , , , 

\K{s) - r"|3 - J 



R(5) 



R(s) - r'|3 



(33) 



with r' a dummy variable. Equation (32) (or equivalently eqs. [33] and [22]) provides the correct expression 
for the force Fa on the satellite moving in the gravitational field of the primary galaxy. 



4.2 Self-gravity of the response 



The shift of the stellar barycenter was viewed by Hernquist and Weinberg (1989) as a manifestation 
of the self-gravity of the response. According to our analysis, the modified kernel in equation (32) derived 
to account for this displacement is introduced to restore the correlation cancelled in the formalism of TLR 
when /°P is introduced to represent the equilibrium stellar system. This is a key correction that accounts 
for the self-gravity of the stellar background. We underline the fact that in TLR the distribution function is 
perturbed only by the satellite external potential. The self gravity of the response has a different origin (White 
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1983; Nelson & Trcmaino 1997): It expresses the "change" in the galaxy's self-interaction potential induced 
by the perturbed density field (excited by the intruder). This would result in an additional component to 
the force Fa on the satellite which is of higher order, and for this reason is neglected in our scheme. In fact, 
the response function A/ modified to account for the perturbed self-gravity would be of order , leading 
to correction of order on the drag force. 

This considerations are consistent with results, from numerical simulations by Bontekoe & Van Albada 
(1987) and Zaritsky & White (1988), on the negligible role played by the perturbed self-gravity in affecting 
the process of energy and angular momentum transfer by the satellite. The self-gravity of the equilibrium 
system is included and affects the process of energy transfer (as it will be shown in §6.5). 



V. FRICTION AND TIDAL FORCE DECOMPOSITION 



5.1 High speed encounters 



In this Section we examine the case of a satellite M moving at high speed relative to a spherically 
symmetric stellar background. In this limiting case and to order "G^", the dynamics of the satellite and 
that of the stars can be approximated as 

r(s) = r - (i - s)v (34) 

R(s) = R - (t - s)V, (35) 

with V and v constant vectors. The satellite and the stars thus move along straight lines, and we neglect the 
self-gravity of the stellar equilibrium system. The validity of this approximation thus restricts to encounters 
whose duration docs not exceed the crossing time of the stars in the galaxy. In the high speed limit, the 
density disturbance thus reads 



An(r) =GMm(3 f ^ 
\ 2tt 



3/2 



/ ds f dsv no(|r - {t - s)v|)e[- 

J —OO J 



R - (t - s)V -v+{t-s)v 
|R - (i - s)V - r (t - s)v|3 



R - (t - s) V - r' 



s)V-r'|3_ 



(36) 



as the presence of the Dirac function 53(r"(< — s) — r) in the coordinate space relates the density at r to the 
density at [r — (i — s)v]. The spherical symmetry is lost since the satellite produces a distortion that keeps 
memory of the star and satellite motions; as a consequence, An is a function of vector r. 

The density response, when V ^ a, should not depend on a explicitly and its expression can be estimated 
in the limit [3 ^ oo. The relevant contributions to An arise, in this limit, from those stars having very low 
velocities for which (jv'^ is finite. Hence, to 0{1/V'^), the density response An can be evaluated expanding 
the arguments of the integrals to first order in the velocity v (the calculation is sketched in Appendix B). 

We find that the response function An(r) naturally decomposes into two contributions: the first repre- 
senting the wake rised by the satellite along his trail: The stream of particles behind M is characterized by 
an overdensity 

/■oo 

Ain(r) = A-KGMno{r) I dr t^3(R - r - tV). (37) 
Jo 

The second contribution rises from the density inhomogeneities of the underlying unperturbed stellar system. 
Introducing the unit vector n = Y /V, we have 

A2n(r) = -^D(R_r) • Vrno(r) (38) 
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with 



Dy = S + n 



'^^ (|Y|)-Y.n 
R-n-n 



Y - n(Y • n) _ R - n(R n) 
^ (|Y|)-Y-n R-R n ' ^ > 



and 



The distortion vector D gives the description of the global tides excited by the satellite during the encounter. 
In equation (39) the vector S is defined as 

S(R)= / dxx-^ ^n(|R-nx|) (40) 

Jo |R ~ '^•^l 

n{r)= f d^r'noir') (41) 

J r'>r 

We have verified that the center of mass of the perturbed stellar system is not affected by the density 
distortion An = Ain + A2n. The displacement vector S compensates for the shift of the center of mass 
caused by the global tide, in the frame comoving with the galaxy. The result illustrates the importance 
of defining the appropriate physical frame of reference for the description of the gravitational interaction 
between the satellite and the primary stellar system. 

Given the density distortion An, we can calculate the force acting on the satellite. We find that Fa 
results in the superposition of two components., Fa = F^p + Ft, the first representing the frictional drag 
induced by the density enhancement Ain that originates behind M 

F.. = -n4.[^^J^ r.,P0(|R-rV 



r ..M^:^ (42) 

Je/V T 



where the stellar mass density is defined as pQ{r) = Nmno{r). In the expression of the force, the Coulomb 
logarithm (eq. [25]) is replaced by an integral, i.e., an "effective" In A which depends on the density distri- 
bution along the wake, since the vector R — rV selects only those stars streaming behind the satellite in its 
motion. We thus find that dynamical friction rises from a memory effect involving only those stars perturbed 
along the motion of the satellite. We may consider, loosely speaking, this contribution to the force as local 
since the drag on M is induced by the bow shock that the satellite excites in its supersonic motion across 
the medium (see eq. [26]). Equation (42) reproduces Chandrasekhar formula for an infinite homogeneous 
medium. 

In addition to the frictional force, we discovered a new component that is closely related to the tidal 
distortion induced by the satellite in its passage. To first order in the perturbation, this component is 

[GM]2 /■ , ^ „ , , R 



1/2 



The force Ft involves a volume integral (contrary to Fnp) depending on the density gradients of the 
unperturbed stellar system: It therefore accounts for the nonlocal response of the galaxy to the disturbance 
and can be viewed as a force resulting from the tidal deformation rised by the satellite during the gravitational 
interaction. We will refer to as tidal force thereafter. Because of the nonhomogeneity of the underlying 
stellar background. Ft is a force acting not exclusively along the direction of motion: In a spherical system 
it acquires a component along R, i.e., along the vector joining the center of the galaxy to the instantaneous 
position of the satellite. Notice that the decomposition of the force in terms of a tidal and frictional drag 
has been possible since we have described the interaction of the spherical system from the noninertial frame 
of reference comoving with the center of mass of the galaxy: The galaxy, as viewed in the laboratory frame, 
suffers in addition the recoil. 



VI. DYNAMICAL FRICTION IN FINITE SIZE CLOUDS 



6.1. Uniform clouds 
11 



As a first application, wc estimate the frictional force experienced by a satellite traveling through a ho- 
mogeneous spherical cloud of radius L. Consistently with equation (35), the satellite motion is approximated 
as uniform and introducing the impact parameter b defined so that b • V = we have 



R(t) = b + tV. (44) 

If V lies along the x-axis, and tp denotes the angle between V and R (varying in between (0, tt)), the satellite 
is seen to enter (exit) the stellar distribution at an angle V'l (V'2) for which 

cot^^i = -i(L2- 62)1/2 (45) 

cotV2 = +^(i^'-fe')^/' (46) 

In this case, the logarithmic integral that appears in equation (42) is readily estimated. When the satellite 
moves inside the cloud (i.e., when '^1 < ^ < </^2 ) it experiences a frictional force 

Fbf = -n 47r-^-|^po InAin (47) 

where 

, , , 6c0tV' + (i2_ 52)1/2^ 

In Ain = In '- — . (48) 

In Figure 1, the effective Coulomb logarithm (In A; solid line) is plotted against il> for two selected values of 
h/L and for a cusp e/L = GM/LV^ = 0.01. We notice that the frictional force, and in turn InA increases 
in magnitude as the massive particle progresses along its motion inside the cloud: Early in the encounter 
the force is small but as the wake develops the drag rises. For sufficiently small (large) impact parameters 
relative to L, the force stays constant (steadily increases) until the edge of the cloud is reached. For b <^ L 
the Coulomb logarithm becomes nearly independent of b and bmax ^ L, the size of the cloud. 

At angles tp < V-'2 the satellite finds itself outside of the stellar distribution. Despite that, the drag force 
is non vanishing because the stream of particles excited during its passage still produces a deceleration. This 
is a consequence of the memory effect present in the correlation tensor. The contribution to the force is 
given by 

^ [GM]^ , bcotf+{L- -Ij')'^^ +e 

with ip < V'2- For (>' 0, i.e., for t 00, the logarithm slowly vanishes as {2/b){L'^ — b'^Y^'^ tan-)/) oc 1/t. 

In a real system, the force will decay as soon as the stars return to equilibrium, after a time t ^ l/w, 
where with u) we denote the mean internal frequency (w =< >^^^ / < >^^^). In the high speed limit, 
stars respond to the perturbation induced by the satellite as dust particles. Hence equation (49) is valid as 
long as the stellar dynamical time exceeds the time scale of the encounter Te, or equivalently ujR/V ^ 1. 

Notice that in addition to the frictional force Fdf, TLR predicts also a tidal contribution that we will 
evaluate in §6.3. 



6.2 Nonuniform spherical galaxy 



In a spherical system with unperturbed density profile 



Po{r) = Po 



(50) 
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the satellite approaehes initially regions of increasing density and consequently suffers a drag whose extent 
with time (i.e., with decreasing At ^ < 7r/2 the medium begins to rarefy and the braking force weakens 
as the wake develops in a medium of decreasing density. 

The magnitude of the force as a function of the orbital phase ip can be computed analytically from 
equation (42) for 7 = 3/2, defining the "effective" Coulomb logarithm as 

lnA= 

Je/V T 

In this model, the mass within radius r declines logarithmically and the satellite in its motion remains always 
confined within the star's background. Figure 2 shows In A for two values of the impact parameter b = 0.9 
and 0.5, (respectively) expressed in units of L, the radial scale below which the density profile is almost 
uniform. We find that if the satellite moves in the halo of declining density only grazing the uniform density 
core, the drag reduces ( by a factor of ~ 2 for & = 0.9) relative to the case of a cloud of radius L, when 
evaluated at phase tp 7r/2. Instead, at small impact parameters {b ^ 0.1) the Coulomb logarithm becomes 
comparable to that of a cloud of uniform density pQ and size L (Figure 1, solid line). 

6.3 Tidal Distortion 

Since the vector (GM/y^)D describes the distortion in the density field induced by the satellite during 
its motion we can draw the isodensity contours of the perturbed stellar system. The first series of plots 
describe the deformation of a density contour lying inside R; the second reproduces the deformation induced 
by the passage of the perturber moving within the isodensity contour in the case of a uniform finite size 
cloud. 

Figure 3 illustrates the evolution of a isodensity level in the orbital plane, for b/L = 1.5. The satellite is 
moving along the trajectory given by equation (44) and the displacement vector is computed for GM/bV"^ = 
0.01. The satellite along its trail rises a tide and a bulge forms which lag behind. 

In Figure 4, we draw a sequence of contour levels describing the extent of the deformation when the 
satellite hits the contour itself: here, b/Lc = 0.5 and GM/bV^ = 0.01. The sharp decrement in the density 
behind M represents the region which is adjacent to the overdensity responsible to the frictional drag. 

6.4 Tidal force in shortlived distant encounters 

At high relative speeds {V ^< >^/^) and large impact parameters {b ^< >^/^, the distortion 
induced on the stellar system acts back to cause tidal drag. By expanding equation [43] in multipoles, we 
find, to leading order. 

Fa = Ft = [GM^Nm^^ (^5{K • n) - 2-^) (52) 

resulting in a total energy loss AE which coincides with that of Spitzer, once equation (44) is adopted for 
the motion of the satellite. The force (52) represents the effect of the global tides on the satellite and has a 
component along V and a component along R, due to the nonuniformities excited by the tidal interaction. 
The radial component is not conservative, i.e., it can not be written in terms of an effective potential but 
affects the energy balance equation. 

In a flyby, the extent of the distortion increases with time. The bulge becomes pronounced after the 
satellite has reached the distance of closet approach and attains its maximum at a later time t ~ l/w; 
the stars will then adjust into a new equilibrium. The results thus seem to indicate that there is a tim,6 
delay in between the moment of closet approach and the moment of maximum deformation of the bulge. The 
magnitude of the distortion D can be computed using equation (39) and is approximately {GM /V'^){D /2L) ^ 
iyi < >^^'^){e/b){b/ L — 1)~^/^. where we indicate with 2L the dimension of the galaxy. 



|R-tVI 



(51) 
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6.5 Energy dissipation 



To gain insight into the new result we need to explore further the role played by the tidal force. Can 
the global tidal distortion depicted in Figures (4) and (5) induce a drag on M , i.e., an effective energy and 
angular momentum loss for the satellite ? How does this loss compare in magnitude to the energy loss by 
friction ? 

We thus need to estimate the total energy change (resulting from the friction and tidal components) 
along the satellite motion. The calculation, outlined in Appendix C, shows that AE < for any arbitrary 
stellar distribution, implying that orbital energy is transferred to the stars, resulting in an effective drag for 
the satellite. A conservative term from the tidal field is found which does not contribute to AE. The total 
energy loss for M decomposes into two terms 



AE = AE+ + AE:,. 



(53) 



The first reads 



A^4 



^ ' -Nm / dt 



V 



B?{t) Jv 



d^r no{r)- 



R(t) • r 
R(i) - r|' 



(54) 



where the domain V is limited to those stellar encounters having |R(t) — r| > e (we here restore the notation 
indicating the dependence of R on time t for seek of clearness). The second is 



A^* = -[GM]'^Nm— 



+ 00 



dt 



+ 00 



dt 



R^t) 



(55) 



where Sit is a function of the actual position of the satellite 



r'<R{t) 



dsv'noir'); 



(56) 



Clt is proportional to the stellar mass comprised within a sphere of radius R{t), and SI = 1 — SI. AE^, gives 
a nonvanishing contribution to the energy loss only if the satellite in its past traveled inside the stellar 
background. 

Restricting the analysis to the case of a finite size cloud of radius L, we are able to show that for b/L > 1 
tidal drag leads to an energy loss given by equation (54) 



AE ■ 



-Sn—^po 



(57) 



In a close encounter (i.e., with b < L) the frictional force becomes important and the total energy loss 
can be written as sum of two contributions: 



AE+ = - 47ri^po 



62 



M 



3Lm In + 3Lm In Lm / e 



and 



1 - 



Lm 
8L3 



(58) 



(59) 



The last term in the square bracket of equation (58) represents the energy loss by dynamical friction 
(from eq. [47] and [49]): here we show that the relevant scale in the Coulomb logarithm is 



bmax = Lm = 2(L^ — 6^)^/^, 



(60) 
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the effective length traveled by satellite during its high speed encoimter, within the galaxy. This is the 
first analytic self-consistent derivation of the maximum impact parameter bmax- The derivation of equations 
(57-59) is outlined in Appendix C. 

Figure 5 depicts the relative contributions to the total energy loss. When b > L the only energy loss is 
induces by the tidal response and is given by equation (57) (solid line). For b < L, the dot-dashed line refers 
to the contribution resulting solely from "friction" (from eq.[42]). The dashed line denotes the contribution 
from equation (59), while solid line gives the total energy loss. Clearly, in equation (58) the additional terms 
resulting from the tidal response tend to reduce the contribution resulting from friction. This goes in the 
direction indicated by numerical findings of Hernquist & Weinberg (1989; despite the misuse of the definition 
of the self-gravity of the response) and Zaritsky & White (1988). The decrement increases with decreasing 
b/L, as illustrated in Figure 5. 

CONCLUSIONS 

We have derived a formalism, within the linear response theory, for the analysis of the interaction of 
a massive particle (the perturber) with a spherical galaxy whose equilibrium is described by a one-particle 
distribution function. 

We applied the pcrturbative technique to the case of a satellite moving at high speed across a stellar sys- 
tem. We found a natural decomposition of the back-reaction force into two components: a global component 
resulting from the "tidal" interaction and a component that is related to dynamical friction. 

The results are relevant only to the description of shortlived encounters since, in the application, we 
neglected the self-gravity of the stellar background. In paper II (Colpi 1997) we will explore the accretion 
of a satellite orbiting around a spherical companion galaxy. In a consistent manner we will account for the 
self-gravity of stars and for the actual dynamics of the satellite. We will present a semianalytical method for 
the description of the process of tidal capture of a satellite, and for the analysis of the evolution of a " binary" 
with the aim at exploring the nature of the gravitational interaction between the two systems. Irreversible 
processes can cause orbital decay by exciting stars in "near" resonance with the orbit of the satellite; this 
phase precedes the final plunge in. The drag experienced by a satellite moving along an arbitrary orbit inside 
the galaxy will be the nc^xt problem that we wish to consider. 

During final preparation of this work. Nelson & Tremaine (1997) submitted a paper where a formal 
study of the linear response theory is presented. The authors derive properties of the response operator in 
inhomogeneous dynamical systems and express the cociffieient of dynamical friction in terms of the equilibrium 
fluctuations, analogously to our equation (16). Their theory is developed in an erbitrary inertial frame. We 
have instead chosen to describe tha process of irreversible energy transfer between the two components 
in a frame comoving with the center of mass of the galaxy (eq.[32]). As previously discussed, we expect 
that our formalism is preferable when an effective one-particle dynamics is adopted for the description of 
the unperturbed stellar motion because it automatically takes into account for the conservation of linear 
momentum in the scattering process of the galaxy whit the satellite. 

We gratefully thank B. Bertotti, F. Combes, G. Gyuk, and I. Wasserman for useful discussions. This 
research was carried out with financial support from the Italian Ministero deU'Universita e della Ricerca 
Scientifica e Tecnologica. 



APPENDIX A 



1. Total energy loss 

In this Appendix we sketch the derivation of equation (23) for the energy loss experienced by the satellite 
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during the encounter, using equation (19): 

/ + 00 l-t !• 

dt / ds I dzTdsw f°'^vlT'"'Vt^ 
-OO J —OO J 



{Al) 



where index t (or s) denote that the variable is evaluated at time t (or s). The domain of integration over 
time s and t can be interchanged to give, in an equivalent form 



/ + 00 P + OO r- 

ds dt d^rdsv r^v^ T^" {V + v°-- v") . 

-OO J s J 

Restoring the complete expression for the correlation tensor we have 
AE = - [GMfNrn^li j d^r d^w /°p 



{A2) 



ds 



dtvl 



R..S 



d 



dt IRt - Ft 



|Rt-rt|3 



(A3) 



In the second term of equation (A3) the vectors R, r and v are evaluated at time t. 

Because of the isotropy of the unperturbed stellar system, the contribution to AS resulting from the 
first term in squared bracket. 



r+°° T}b _ b 



J dsvds^r j 



|R«-r,|4 

vanishes identically. The second term gives instead equation (23) together with (24). 



(A4) 



2. Impulse approximation 



In the impulse approximation, the motion of the satellite and of the background stars is uniform. We 
can evaluate Ai? according to the simplifying assumption that Rj = b + tW and = r + tv, with V and v 
constant vectors and b • V = 0. Defining a = b — r and w = V — v, the scalar J is found to be 

J= V ^ + % =2vB {A5) 

with 

B = 4^. (A6) 

In equation (A6) a_L is the component of a in the orthogonal direction of vector w and is equal to a_L = 
a — (a • w)w/w^. The total energy loss thus reads 

AE = -2[GMfNm^P j dsr ci3v/°P(v • B)^. (AT) 

The energy loss by recoil is recovered from equation (A7) in the limit b oo and V ^ oo valid if the 
satellite moving at high speed maintains always at a distance from the stellar distribution. To first leading 
order, the vector B reduces to B = b/(V^&^) yielding 

AE = -2[GM]'^Nm/{V'^ b^). (AS) 

(We recall that / d3v/°P(w«)2 = a'^no{r).) 
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Instead, in the limit F — > oo but rjh finite, 



B 



a+(r- V)V/F^ 
y[a+(r- V)V/F2]' 



and the energy loss equals 



V)2/y2]- 

Expanding equation (AlO) to order < > /b'^ we recover Spitzer result (eq. [28]). 



{A9) 



(AW) 



APPENDIX B 



Wake and global tides 



In this Appendix we summarize the main steps that lead to the decomposition of the force Fa in terms 
of FnF and Fa- Both components are derived to order 0{1/V'^). 

The density response An is calculated from equation (36), in the high speed limit, according to equations 
(34) and (35). The force term (the one in squares brackets) and the density field no{\r—{t—s)v\) are expanded, 
in equation (36), to first order in the velocity v, contributing to order GM/V"^. 

Let us denote with Ain the distortion resulting from the expansion of the force term, and with 
the one from the expansion of the unperturbed density no- The global distortion will therefore be 



An = Ain + A2n 

For ease of analysis, we introduce the vector 

X = R-r- (<-s)V 

so that 

= ^ + {t-s)v'—^—^+0{v^) 



|X+(t-s)v|3 X3 



(Bl) 

m 
m 



to first order in v (X has no relation with the one introduced in §4.1). 
According to equation (36) 



Am =GMm(3 ' j ds J dgv no{r)e-i'^"'''' 

'X" R''-{t-s)V'' 



X3 |R_(f_s)V|3 



d^r' no(r') 



r'<d 



+ {t-s) V'' 



d 

dx^lc^ 



where d simply denotes |R — (t — s)V|. The term in squared braked vanishes identically because of isotropy. 
Only the term oc v"" gives a contribution to the distortion. Since for a Maxwellian distribution function 



/ 



rfsv e" 



r3/2 



-m(3 



-5/2 



io6 



we finally have 



Ain(r) = GMn, 



o(r) /* 

J — C 



ds {t - s)Vx 



|X|3 
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Using the divergence theorem we have 

Ain(r) = 47r GMno{r) f ds {t - s)(53(X) 

J — OO 



m 



which is equivalent to equation (37). 

The global distortion in the density field is derived expanding no to first order in v; using Gauss theorem 
we derive the following expression 



A2n(r) =GMmP ( ^ ) ( 



m[3 



3/2 



dnr 



v" v'' 



2tt 



ds (t — s) J dave ^ 
dsr' no(r') 



{B8) 



X3 |R-(t-s)V|3^^,<, 

Using equation (B5) and introducing the velocity vector n = V/V, we find that the density response reads 

GM dn 



A2n(r) 



r+oo 



y2 Q^a 



dx X 



R-r-nx^ R-na;3 R-r-na;^ ^' 



(B9) 



where fl is defined by equation (41). Equation provides the expression of the displacement vector in integral 
form and after lengthy calculations one recovers equation (39) for D. 



APPENDIX C 



l.The total energy loss 



In this Appendix we estimate the extent of the energy loss resulting from friction and from the tides 
excited by the satellite during its passage: 



/ + 00 
dtY -Fa (CI) 
-OO 

Within TLR, the macroscopic force on M is found to be 

Fa = -GMiVm J rf3r[Ain(r) + A2n(r)] |^ "^'jg (C2) 
with Ain and A2n given by equations (37) and (38). The contribution coming from Ain directly gives: 

^ _MGMrNm r- ,^no(|R-V.|) ^^3^ 

y Jo 

while the term containing A2n in eq. (C2) can be written as the sum of three terms: 

J '^sr A2n(r) ^— ^ jd,v (C4) 



DH-rno(r)^-^ 
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where Dr-t is given by equation (39). The first term vanishes due to the divergence theorem, the sec- 
ond exactly cancels the Chandrasekhar term (C3) leaving only the third contribution which contains both 
frictional and tidal effects to give: 



{C5) 



The force (C5) can be decomposed into the sum of a conservative contribution plus a component leading to 
dissipation of energy: 



Fa =Vr 



^ ^ Nm J rf3rno(r)D(R_r) ^ ' 



1/2 



Nm / d3rno{r) 



|R-r|3 



|R-r|3 

VRl?^R_r)- 



For ease of analysis we introduce the vector 



C = R-r 



(G6) 



expressing the relative distance between the satellite and the stars of the background. In the dissipative 
component the spatial integration is then limited to a domain V satisfying the condition ^ > e (see eq. [27]). 
The vector D 



n*'ln(e-e-n) + 



f^-n'X^-n) 



can be schematically written in the form 



n''ln(i?-R-n) + 



i?^-n''(R-n) 



i?- R- n 



+ 5" 



where S is defined by equation (40). 

In order to calculate the total energy loss we need to evaluate explicitly the component of the force 
along the velocity unit vector n. Along the unit vector n, the force is given 



Nm j dsr no{r)n' 



One can prove that 



and 



[GM]' 



^3 dR" 



^3 QJla L"' ^3 ^(^ _ n • R) 

Substituting (Cll) and (C12) in (CIO), and using Gauss theorem, we find 



Noting that 



n{R) 



[GM]' 



■Nm 



n(R) 



nojr) 
R-r|3' 

R r 



± J dsr no(r) = J d^r^^^ ' ^2 J ^sr no(r)^ 



i?3 ^- 

we find that AE splits into two terms AE+ + AE* given by 

AE+ = .J Nm 



V 



rfsr no(r)- 



R r 

R-r|3 



(CIO) 

(Gil) 
(G12) 

(CIS) 
(CM) 

(G15) 
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and 



V 



(C16) 



Both terms are negative as shown below. 

As regard to the first term AE+ (C15), owing to the isotropy of the background, we can arbitrarily 
select a direction for R (along z axis) so that the resulting volume integral 



/ 



dsT no(r) 



Rz 



(i?2 +r2)3/2 



[l-2ii^/(i?2 + r2)]-3/2 



(C17) 



has a kernel of the form (1 — z)~^^^ = [1 + 3z/2 — 15z'^/8 + 0{z^)] when expanded in series. Since only 
the terms with odd powers give a nonvanishing contribution, the volume integral has definite positive sign, 
resulting in a net energy loss. Also the term in equation (CI 6) can be suitably written in a form involving 
only positive functions. 

We know evaluate AE*. If we explicitly insert the expression of S as given by equations (40) and (41) 
into equation (C16) we have that the integral 



is equal to 
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giving a contribution to the energy loss 



A^, = -[GM]^Nm dt^^^ 
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where we recovered the dependence upon t ofHt, denoting the satellite position at current time t. 

Since, in a high speed encounter, the satellite moves along a trajectory that can be approximated as 
linear, we can evaluate explicitly the integral over the variable s using equation (44). Noting that the 
following equality holds 
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after some calculation we find that 
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Using equation (C21) we recover eq.(55). 



2. Energy loss in a uniform cloud 



Here we provide a quantitative estimate of the total energy loss and determine the relative importance 
of the tidal and frictional contributions to the process of dissipation. We consider the case of a uniform cloud 

of mass density po ~ Nmno and radius L. 

CASE A: If the satellite travels outside the spherical cloud the only contribution to the energy loss comes 
from the tidal field and is given by 

AE+ = -^^Po / dt 

' J — oo 



3 -Rj Rt Rt 



(G23) 



20 



Since Rt = (6^ + (Ft)^)^/^, the integral over time t can be evaluated analytically yielding 



AE+ = -Sn^^po 



L-g-(6^-XY/^sin-f- 



(C24) 



where b is the impact parameter {b> L). In the limit L equation (C24) reduces to (28). 

CASE B: If during the shortlived encounter the satellite enters the spherical cloud (6 < L), the total energy 
loss is given by equation (54) plus (55). Notice that the satellite travels a distance Lm = 2(L^ — 6^)-^/^, 
within the cloud, over a time = Lm/V . In this case 



[GMf 
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(C25) 



where 7< (7>) is the volume integral resulting from the region of the cloud with [0 < r < R — e] 
{[L > r > R + e] respectively): 

167r „ , L^-R? 
/<+/>= — + 27rln 



while 



47rL3 L ^ R + L 



The time integrals can then be carried out analytically: 



(C26) 
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Similarity, we can carry out the integration of equation (55) involving the functions O and CI. The 
calculation is simple and lead to the following expression for 



4^ [GMf Ll^ 
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where Lm = 1\fl?^-W is the length of the portion of trajectory within the cloud. 
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Figure Caption 



Figure 1: Effective Coulomb logarithm versus phase ?/; for a uniform finite size cloud Solidline is for b = 0.1; 
dashed line for b — 0.5 and dashed-dotted line for b — 0.9. 

Figure 2: Effective Coulomb logarithm {solid line; eq.[51]) against phase ip ioi b — 0.9 (right panel) and for 
b = 0.5 (left panel); e/L — 0.01. The spherical stellar background has density profile given by eq. [50] with 
7 — 3/2. Dashed line gives the comparison with the case of a uniform density cloud. 

Figure 3 Isodensity contours, in the orbital plane {x,y) for b/L — 1.5 : Square denotes the position of the 
satellite moving along the x axis, from left to right. The phase ip decreases as time elapses and the four 
panels refer to ip — 120, 90, 60, 45 degrees respectively. 

Figure 4 Isodensity contours, in the orbital plane {x, y) for b/L — 0.5 : Square denotes the position of the 
satellite moving along the x axis, from left to right. The phase ^ decreases as time elapses and the four 
panels refer to V' = 157, 130, 50, 26 degrees respectively. 

Figure 5 Energy loss — Ai? in units of [GMj^po L/V"^ as a function of b/L. The spherical system is uniform 
and has size L. For b/L > 1 only tides cause dissipation of energy (eq. [57]-sohd line) For b/L < 1 the tidal 
contribution of equation (59) is indicated by the dashed and the dot-dashed line denotes the contribution 
from the logarithmic term of equation (58) with e/L = 0.01. The solid line gives the total energy loss. 
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